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Abstract. We are familiar with Dirac equation in flat space by which we can 
investigate the behaviour of half-integral spin particle. With the introduction of 
general relativistic effects the form of the Dirac equation will be modified. For 
^ ' the cases of different background geometry like Kerr, Schwarzschild etc. the cor- 

responding form of the Dirac equation as well as the solution will be different. 
f— ^ \ In 1972, Teukolsky wrote the Dirac equation in Kerr geometry. Chandrasekhar 

^^ ' separated it into radial and angular parts in 1976. Later Chakrabarti solved the 

angular equation in 1984. In 1999 Mukhopadhyay and Chakrabarti have solved 
the radial Dirac equation in Kerr geometry in a spatially complete manner. In 
this review we will discuss these developments systematically and present some 
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1. Introduction 

Behaviour of particles with half integral spin can be investigated through 
the study of Dirac equation. Generally, we are familiar with the Dirac equa- 
tion and its solution when the space-time is flat. In the curved space-time 
where the influence of the gravity is introduced, the corresponding equation 
will be changed in form. Its solution will also be different. In 1972, Teukol- 
sky [1] wrote the Dirac equation in curved space-time particularly in Kerr 
geometry [2] using Newman- Penrose formalism [3]. Through this modified 
Dirac equation we can study the behaviour of spin half particles around the 
spinning black holes. Due to presence of central black hole the space-time 
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is influenced and behaviour of tlie particle is changed with respect to that 
of flat space. ^From the same equation of Teukolsky, Dirac equation for 
Schwarzschild metric [2] (Schwarzschild geometry), where the central black 
hole is static can be studied just by putting the angular momentum pa- 
rameter a of the black hole to zero. So one can study how the behaviour of 
spin half particle in curved space time is influenced by the angular momen- 
tum of black hole. In 1976, Chandrasekhar [4] separated the Dirac equation 
in Kerr geometry into radial and angular parts and solved the radial part 
of the equation asymptotically. Chakrabarti in 1984 [5] solved the angular 
part analytically. Here we shall introduce the spatially complete analyt- 
ical solution of radial Dirac equation [6-7]. So the the complete solution 
of Dirac equation can be studied. Far away from the black hole the the 
modified Dirac equation for curved space-time (for Kerr and Schwarzschild 
geometry [2-3]) and its solution reduce into that of the flat space. 

In this review we will first indicate how Dirac equation in curved space- 
time can be written using Newman-Penrose formalism [3]. Newman- Penrose 
formalism is one of the tetrad formalism where null basis are chosen instead 
of orthonormal basis. To fulfill the understanding of Dirac equation in this 
formalism we also need to know the 'Spinor Analysis' [3]. In the next Sec- 
tion, we will briefly describe this in the context of our present purpose. In 
§3 we will write the Dirac equation in Newman-Penrose formalism for flat 
and curved space-time. For curved space we will separate the Dirac equa- 
tion under the background of Kerr geometry. In §4 and §5 we will briefly 
outline the angular and radial solution of Dirac equation respectively. In §6 
we make concluding remarks. 

2. Spinor Analysis 

In Minkowski space we consider a point x* (i = 0,1,2,3) on a null ray 
whose norm is defined as 

Now, we consider two complex numbers ^° and ^^, and their complex con- 
jugate numbers ^^ and f^ in terms of which each point can be written 
as, 

xi = i=(^°r+e^r) (26) 

x' = ^{ee'-ee') (2c) 
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x' = -^{ef-ee') m 



Conversely, we can write, 

1 

1 

1 
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^0^-0'= 1 (2,0^^3) (3^) 



^0^1'^ 1(^1^ .^2^ (36) 



^1^-0'= 1 (^l_,^2) ^3^) 



e^f = ^(x°-x3) (3d) 



Let, 

C^ = a^B^^, (4a) 

e?' = «^:e^' (46) 

where, {A, B, A' , B' = 0, 1), are the hnear transformations in complex two- 
dimensional spaces. The transformation of x* is defined as, 

xl = ^ixK (5) 

Now, using equation (2) and (3) we can write, 

xl = ^{aie + oiie){4'~e' + 4e') + ^ne + a\e){4e' + «};r ) 

= 2("o"o' + ao"o')(a;° + x^) + -(a? a?/ + a\aX>){x^ - x^) 

+-yoiQaii + aQaii){x +ix ) + -(a^QQ/ + aiao')(^ ~ ^x ). (6) 

Similarly, we can write down the relations between x\, xl and x\ with a's 
and x's. Therefore, keeping in mind (5) we can write. 
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Now, imposing the condition that the transformation (5) is Lorentzian we 
can write, 



1 



So, 



This gives, 



0— n' 1— V 0—0' 1—1' 

a^a^, + aoQQ/ aQa]^, + aQai, 

0— n' 1—1' D— n' 1—1' 



AA = 1 



1. (7) 



(8) 

Now we consider A = A = 1, so individually each transformation of ^ 
is Lorentzian. So we can conclude if transformation (5) is Lorentzian, the 
necessary condition is transformation (4) is also Lorentzian. 

Now we define spinors ^ , t] of rank one as £,^ = O-q^^ and 7]^ = 

1. Since ^"^ and r]^ are 



a 



4',ry^', {A, A', B, B' = 0), where 



two spinors of same class. 



< 


= 


a% 






tO^l 



1^0 



cw - ert 



which is invariant under unimodular transformation, i.e., 

^AB^ 1] -^ invariant 



(9) 



(10) 



where, cab is Levi-Civita symbol. Here as in the case of tensor analysis cab 
and eA'B' are used to lower the spinor indices as, S,a = ^'^'^CA- 

Now, using above information the representation of position vector x* 
can be written as 



^1^-0' ^1^-1' 



1 
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n "^ 

X + X 

1 ■ 2 

X — IX 



1 • 9 

X + IX 



(11) 



Generally any vector X* can be written in terms of spinor of rank two as, 



X' 



too' 

^10' 



01' 

11' 



1 
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X" + X^ X^ + iX^ 
X^ - iX"^ X^ - X^ 



X 



AB' 



(12) 



So a 4- vector is associated with a hermitian matrix such that, 
.^-,2 .(xi)2-(x2)2-(x3)2 = {X^+X^){X^-X^)-{X^+iX'^){X^-iX'^) 



{X^ 



^00' til' tOl', 



tOO'( 



^11' 



-10' I c01'( 



2(r e'^ -e'^'^io') = {rim'+iive'+iwe'+e'i.v 

of norms, we can write it in 
gijX^X^ = eAC^B'D'X ^ 



Xab'X 



AB' 



^From the definition of norms, we can write it in two different representa- 
tions: 

-AB' vCD' 



(13) 
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Therefore, we can transform X* ^^ X 



using, 
AB' 



(14a) 



X 



AB' 



AB' Y'i 



(146) 



where, a/^ and cr\^i are nothing but Pauli matrices and their conjugate 
matrices with a factor —7=. 

A curved space-time is locahy Minkowskian. At each point of space-time 
an orthonormal Dyad basis can be set up as C,f\ and C,^,^ (a, a' = 0, 1 and 
A, A' = 0, 1) for spinors. 

We define, Cfg-, = o and Cn~) = ' • The condition of orthonormahty is 



eABO^l^ = o'^l^ - oH'' = obI^ 



-oHa = 1. 



Also it is clear that, e^^^^^^^l^Cf) = e 



AB 



/^F. 



Then the null vectors are defined as 

Where, vectors obey relations of null tetrads such as, 
l^rii = l,m^fhi = —1 and all other products give zero. 
In this way using dyad basis we can set up four null vectors which are 

basis of Newman-Penrose formalism. Using (14) we can write the basis 

explicitly as 

^ = CrAB'C(0)C(0') = (^AB'O O , (15a) 



m 



i /-A aB' _ i ^AJB' 
'^AB''^{0)'^{1') — '^AB'O ' , 



^i _ ^i /-A 7B' _ i ,A-B' 
^ — '^AB'^1)^{0') — ^AB''' O , 

„i — ^i M 7B' _ i jATB' 



(156) 
(15c) 
(15d) 



Thus, in Newman-Penrose formalism the Pauli matrices change their forms 
as. 



a. 



(^AB' 



AB' 



1 

~ V2 


r 


re 


1 


rii 


—fhi 


V2 


- 


-mi 


k 



(16a) 
(166) 



Therefore in this basis, the directional derivatives can be written as, 
D = rdi, A = n*9i, 5 = m^di and 5* = frfdi. 
Thus, the spinor equivalents of Newman-Penrose formalism are 
aoo' = -D, div = A, 5oi' = (5, 9io' = 5*. 
Due to the reason, as explained earlier Vj ^^ ^ ab'-, so we can write. 



ViX, 



X 



j;* 



^ AB'^CD' 



X, 



CD':AB' 



6 BANIBRATA MUKHOPADHYAY 

therefore, 

XcD';AB' = '^CD''^AB'^j;i- (17) 

For covariant derivatives spin coefficients T are introduced. In the Newman- 
Penrose formahsm these different coefficients are assigned in terms of special 
symbols which are given in [3]. 

3. Dirac Equation and its Separation 

Before going into discussion, we should mention about the unit of the sys- 
tem. Here we have chosen throughout h = c = G = 1, where h = Plank 
constant, c = speed of light and G = gravitational constant. It is very clear 
that simultaneously all these quantities are chosen as unity implying the 
corresponding system is dimensionless. 

The Dirac equation in flat space using Newman-Penrose formalism can 
be written as, 

a'AB'diP"^ + ifi*QB' = (18a) 

<^AB'diQ^ + ili*PB' = 0. (186) 

Here, P and Q are the pair of spinors. ^*/v2 is the mass of the particles 
and (t\q, is nothing but Pauli matrix, because l/\/2 factors are canceled 
in the equation. 

In curved space time Dirac equation reduces to 

4B'^;f + i^^*Q^'^C'B' = 0, (19a) 

cT^B'Qti + i/i*^^'ec'B' = 0, (195) 

where, o"^^/ is same as defined in equation (16a). 
Now, consider B' = 0, then (19a) reduces to 

or, 

(aoo'P° + r^oo'^') + {dwP^ + rHio'^') - i^^*Q^' = o, 

Therefore, 

{D + Piooo' - Tooio')^" + {5* + Pnoo' - Pono')^' - i/^*Q'' = (20) 

Similarly, choosing S' = 1, we can get another similar type equation and 
then we can get corresponding conjugate equation of both by interchanging 
P and Q. Now choosing, 

Fi = po, F2 = P\ Gi = Qi', G2 = -0°' 



/ sVp'a 










2aMr/p^A \ 





-A/p^ 
















-l/p' 







V 2aMr/p'^A 








-(A- 


- a^sin^e)/p^Asiv?e J 
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and replacing various spin coefficients by their named symbols [3] we 
get the the Dirac equation in Newman-Penrose formalism in its reduced 
form as, 

{D + e- p)Fi + {5*+7r- a)F2 = i/i*Gi, (21a) 

(A + ^ - 7)F2 + (5 + /3 - t)Fi = ip,G2, (216) 

{D + e* - p*)G2 -{S + TT*- a*)Gi = i/i^Fs, (21c) 

(A + M* - 7*)Gi - {6* +(3*- T*)G2 = ip,Fi. {21d) 



3.1. BASIS VECTORS OF NEWMAN-PENROSE FORMALISM IN TERMS 
OF KERR GEOMETRY 

The contravariant form of Kerr metric is given as [3] , 



9'' 



(22) 
where, E is the energy, a is specific angular momentum of the black hole, 
M = mass of the black hole, p"^ = r'^ + a^cos^O (should not confuse with 
the spin coefficient r(o)(o)(i)(o') = P)i ^^ = (^^ + o^)^ ~ a^Asin^O, A = 
r"^ + a? - 2Mr. 

In Kerr geometry, the tangent vectors of null geodesies are: ^ = ^^ ^ ' E, 
-T- = ibi?, ^ = 0, ^ = -J-E, where r is the proper time (not to be confused 
with spin coefficient r(o){o)(i){i') = ''")• 

Now, the basis of Newman-Penrose formalism can be defined in Kerr 
geometry as (in tetrad form), 

li = —(A, -p^ 0, -aAsin^6l), (23a) 

m = -i^(A, p2, 0, -aAsin^e), (236) 

2p^ 

rrii = {iasm9, 0, —p ,—i{r +a )sin^), (23c) 

PV2 

e = ^{r^ + a^A,0,a), (23d) 

n' = ^{r^ + a^-A,0,a), (23e) 
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rrii = _ ^ {iasm9,0,l,icosec6), (23/) 

fhi and m* are nothing but complex conjugates of rrii and m* respectively. 

3.2. SEPARATION OF DIRAC EQUATION INTO RADIAL AND ANGULAR 
PARTS 

It is clear that the basis vectors basically become derivative operators when 
these are applied as tangent vectors to the function eH'^*+»"</') . Here, a 
is the frequency of the particle (not to be confused with spin coefficient 
r(o){o)(o)(i') = '^) ^-iid m is the azimuthal quantum number [3]. 

Therefore, we can write, 

r=D = Vo,n = A = -^Vl m = 5 = ^4, ffi = 5* = ^A- 

where, 

^ r. iK r - M , , 

Vn = dr + — + 2n^^, (24a) 

Pl = 8,.-f + 2,.^, (246) 

i^n = de + Q + ncotO (25a) 

Cl = dg-Q + ncote. (256) 

K = {r'^ + a^)(T + am, Q = aasinO + mcosecO. 

The spin coefficients can be written as combination of basis vectors in 
Newman-Penrose formalism which are now expressed in terms of elements 
of different components of Kerr metric. So we are combining those different 
components of basis vectors in a suitable manner and get the spin coeffi- 
cients as, 

K = a = X = Ty = e = 0. (26a) 



1 o cotd iasind 

iasinO 



^ iasind ,, A ., ,,1 r—M _ _ a* /'o«■k^ 

^ = ~7V2^'^ = ~2^^'T = ^+^^'" = ^-^ • (26b) 
Using the above definitions and results and choosing /i = p*Fi, §2 = 
pG2, f2 = F2, 91 = Gi the Dirac equation is reduced to 

Vofi + 2~^/^£i/2/2 = {ifJ'*r + a^^cosO)gi, (27a) 

APJ/2/2 - 2^l'^C\iJi = -2(i/i,r + aii,cose)g2, (276) 

^052 - 2^^/'^C\i^gi = {ifi^r - afi^cos6)f2, (27c) 

AP}/251 + 2^/^/:i/252 = -2(i/i,r - ati^cos9)fi, {27d) 
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Now we will separate the Dirac equation into radial and angular parts by 

choosing, 

h{r,e) = R^i/2{r)S^i/2i0), f2{r,e) = Ri,2{r)S„2{e), 
gi{r,e) = Ri/2{r)S_y2{0), g2{r,9) = ii-i/sW^i/sW- 
Replacing these fi and gi {i = 1,2) into (27) and using separation 

constant A we get, 

-Ci'S't 1 = —(A — am„cos9)S_i (28a) 

2 '2 2 

£ts'_i =+(A + am„ cos 61)5', 1 (286) 

- 2 '2 



2 

1 _ _ , , . X . 1 

1 
'2 



AiVoR^i = (A + impr)A2R^i, (29o) 

A^vIA^R^i = (A - impr)R_i, (296) 

where, rrip is the mass of the particle which is nothing but 2^'^;U^<. Also, 
2^''^R_i/2 is redefined as R-i/2- 

Equations (28) and (29) are the angular and radial Dirac equation re- 
spectively in coupled form with the separation constant A [3]. 

4. Solution of Angular Dirac Equation 

Decoupling equation (28) we obtain the eigenvalue equation for spin-2 par- 
ticles as 

S_i=0. (30) 



L\L\ -I- -L\ -I- (A —a m„cos 0) 

2 2 A -I- anioj cos '^ - '^ 



p ' 



2 



Similarly, one can obtain decoupled equation for spin -|- 2 particles. Here, the 

separation constant A is considered to be the eigenvalue of the equation. 

The exact solutions of this equation for A and 5_i is possible in terms of 

2 
orbital angular momentum quantum number / and the spin of the particle 

s when the parameter pi = — ^ = 1. When the angular momentum of the 

black hole is zero i.e., Schwarzschild case, the equation is reduced in such a 

form that whose solution is nothing but standard spherical harmonics such 

as [8-9], 

S^i/2ie)e'"''^ =^iYUe,<P), (31) 

the eigenvalue i.e., the separation constant can be solved as, 

A2 = (/ + 1/2)2. (^32) 

Similarly, for spin-|-2 particle one can solve <S'_|_i/2 as 

S+i/2{e)e'"''f'=^YUe,<l>), (33) 
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with same eigenvalue A. 

For any non-integral, massless, spin particle the solutions are [8-9] 

S±Me'""^=±sYiU0,<l>), (34) 

A2 = (/ + |s|)(/-|s| + 1). (35) 

In the case of Kerr geometry, when a ^ the equation can be solved by 
perturbative procedure [5] with perturbative parameter aa. The solution 
for pi = TUp/a = 1 and s = ±^ is [5] 



X^ = {1 + -f + aa{p + 2m) + a^a^ 



where, 
and 



y 



^ '""5 '™+2(Z + l) + ac7x 

c. _ ^ a<yy 

2 2 2(i -I- 1) -I- aax 



2{l + 1) + aax 



lYu 



2 



p = F{l,l); x = F{l + 1,1 + 1); y = F{l,l + l) 



(36) 

(37a) 
(376) 

(38) 



F(/i,/2) = [(2/2 + l)(2/i + 1)]^ < hlmQ\hm > [< hl\Wi\ > 

+ (-1)'^-' < hlmOlhm > [< hl^O\h^ > +i-iy^'^piV2 < l^l-^l\l^^ >]] 

(39) 
with < ....|.. > are the usual Clebsh-Gordon coefficients. 

If /?! 7^ 1 then exact solution is not possible. In those cases the analytic 
expression of eigenvalue and angular wave-function are found as infinite 
series not in a compact form as the case pi = 1 . 

^From the general convergence of series expansions one can truncate 
the infinite series upto certain order for particular values of /, s and m. For 
I 



2' 



"2 and m = — 2, up to third order in aa, one obtains [5], 



A' = {l + -)^ + aafi{l,m) + (aa)^ f2il , m) + (aa)"" hil,m), (40) 



'2 2 2 



sinO — I sin sinOcos — 



"^aail + pi) + —{aa)\l - pD 



+-(aa)\l-pl) 



sin 6sin —cos — \- 3sin—cos — 

2 2 2 2 2 



(41) 
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The accuracy of eigenvalues and eigenfunctions decreases as ao" ^ 1. 

5. Solution of Radial Dirac Equation 

In the radial equation independent variable r is extended from to oo. For 
mathematical simplicity we change the independent variable r to r* as 



2Mr+ + am I a ( r 

r* = r-\ log 1 

r I — r_ V r I 



2Mr_ +am/a, ( r , . . 

^log — -1 42 

r I — r_ \ r_ ' 



(for r > r+), here in new r^< co-ordinate system horizon r+ is shifted to — cx) 
[3], so the region is extended from — oo to oo. We also 



unless o" < 



2Mr4 



choose i?_i = i^_i, Aai? 1 = P, 1. Then we are defining 



2„2\ 



(A lb impr) = exp{-^i6)^J{\ + nipr 



and 



^+1 

^2 



■ih.i exp 

"•"2 



ijj^i exp 



1 

— -i tan 

1 

+— i tan 



-1 



Finally choosing, 



^± = V'+i =t V'.i 



and combining the differential equations (29) we get, 



and 



where. 



and 



_d_ 



4- + w]z. 

ar^ 



iaZ^, 



W 



1 _i fm,„r 



A5(A2 +mlr'^fl'^ 



a;2(A2 + mlr'^) + \mp/\/2a ' 

where, cj^ = — . 

Now decoupling equations (43a-b) we get. 



d' 
df^.' 



+ a^]Z± = V±Z±. 



(43a) 
(436) 



(44) 



(45) 
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where, V± = W"^ it -^ and f* is extended from — oo (horizon) to +00. 

The equation (45) is nothing but one dimensional Schrodinger equation 
[10] with potentials V± and the energy of the particle a"^ (since the system 
is dimensionless) in Cartesian co-ordinate system. The equation (45) can be 
solved by WKB approximation method [10-11]. The corresponding solution 
is [6-7], 

Z± = ^ — exp{iu±) ± „ — exp{—iu±) (46) 

\/k± yk± 

where. 



and 



k^ = J{a^-V^), (47) 



u± = k±df^. (48) 



Now we improve the solution by introducing space dependences on coeffi- 
cients A± and B± [6-7] (this is beyond WKB approximation, because WKB 
deals with solutions with constant coefficients). It is seen that far away from 
a black hole, potential varies very slowly. Thus, in those regions one can 
safely write, 

A±- B± = Constant (= c). (49) 

Since the sum of reflection and transmission coefficients must be unity, 

Al+Bl = k±. (50) 

Near the horizon it is seen that potential height reduces to zero so the re- 
flection in that region is almost zero and transmission is almost 100%. This 
is the inner boundary condition. Solving (49) and (50) we get analytical 
expression of space dependent reflection and transmission coefficients far 
away from the black hole which satisfy outer boundary condition. Combin- 
ing the inner and outer boundary conditions, we get analytical expression 
of space dependent coefficients A± and B± which is valid in whole region 
(— cx) to -|-oo). For details see [6-7]. The space dependency of A± and B± 
i.e. the transmission and reflection coefficients arises due to the variation of 
potential with distance. So from the analytical expressions one can easily 
find out at each point what fraction of incoming matter is going inward and 
what other fraction is going outward as a result of the interaction with the 
black hole. These space dependent transmission and reflection coefficients 
are given below [6-7], 



(ci + f) / c , /- -\ , 2fc± 



^* = g = ''7^l'^' + 2 + \/^'*-'^'j + '^l^ f^'"' 



c 
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hi (c2 - f ) / c r- r\ 2k±-c^ , ,, 

i?± = -± = ^ \ ^' C2 - - + V2A;± - c2 + ^- . (516) 



k± h± \ 2 ^ J 4/i± 

Here, a± and 6-|- are defined as 

A± . . 

a± = , (52a) 

which are transmitted and reflected amphtudes of the solution with modi- 
fied WKB method (going beyond WKB method) and 

hi = (ci + 0' + [c2 - 0' + (ci + C2)v/2fcl^+ ^^^±^, (53) 

where, ci and C2 are two constants introduced to satisfy the inner boundary 
condition. The final form of the solution is 

Z± = — exp{iu±) ± — exp{—iu±). (54) 

^Jk± \/k± 

Since the relation between Z± and R,i is known, one can easily calculate 
the radial wave function Ra_i. 

6. Conclusions 

In this review we write the Dirac equation in curved space-time and par- 
ticularly in Kerr geometry. From this, the behaviour of non-integral spin 
particles can be studied in curved space-time. From the form of the equa- 
tion and its solution it is clear that in curved space the particles behave 
in differently than in a flat space-time. The Newman-Penrose formalism is 
used to write the equation where the basis system is null. Dirac equation 
is separated into angular and radial parts. Similar separation can be pos- 
sible on the background of Dyon black hole [12]. The solution of angular 
component of the Dirac equation is first reviewed. The exact solution is 
possible for Hiz = i^ otherwise the solution is approximate [5]. Unlike in 
the case of a Kerr black hole, the solution of the angular equation around 
a Schwarzschild black hole is independent of the azimuthal or meridional 
angles [5, 13,14]. This is expected because of symmetry of the space-time. 
The radial Dirac equation is solved using WKB approximation more 
clearly modified WKB approximation [6-7], where the space dependent 
transmission and reflection coefficients are calculated. Although WKB method 



14 BANIBRATA MUKHOPADHYAY 

is an approximate method, it is improvised in such a way that spatial de- 
pendence of the coefficients of the wave function is obtained. This way we 
ensure that the analytical solution is closer to the exact solution. The re- 
flection and transmission coefficients were found to distinguish strongly the 
solutions of different rest masses and different energies. The solution might 
be of immense use in the study of the spectrum of particles emitted from 
a black hole horizon (Hawking radiation). 
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